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We investigate the structure and the novel emerging features of the mesonic non-singlet spectrum 
of the Minimal Walking Technicolor (MWT) theory. Precision measurements in the nonsinglet 
pseudoscalar and vector channels are compared to the expectations for an IR-conformal field theory 
and a QCD-like theory. Our results favor a scenario in which MWT is (almost) conformal in the 
infrared, while spontaneous chiral symmetry breaking seems less plausible. 
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I. INTRODUCTION 

The idea of a new strong force in our model of Nature 
to explain Electro-Weak Symmetry Breaking (EWSB) 
dynamically was first suggested many years ago 0, 0. 
Borrowing from our intuition of QCD, a new strong sec- 
tor beyond the standard model was proposed in which 
chiral symmetry breaks down at the TeV scale leading to 
EWSB and providing an explanation for the observed 
gauge boson masses. Moreover, the standard model 
fermions acquire their masses via extended technicolor 
interactions. The first models based on these ideas were 
obtained by a naive rescaling of QCD, i.e. they were 
based on a SU(N) gauge theory with a small number of 
fundamental matter fields. Despite the elegance of the 
proposal, it was soon shown that such models are not 
viable candidates: together with the mass of SM par- 
ticles, large Flavor Changing Neutral Currents (FCNC) 
and large values of the Peskin-Takeuchi 0, 0| parameters 
would also be generated. Electroweak precision tests per- 
formed at LEP [5[ put tight experimental constraints on 
FCNC and the oblique parameters, which are incompat- 
ible with such predictions. 

However, the naive rescaling arguments leading to the 
above conclusions can be flawed if the dynamics of the 
new sector is sufficiently different from QCD. In fact the 
intrinsic difficulty of handling strongly interacting mod- 
els has not stopped the theoretical speculations. Walk- 
ing and conformal technicolor theories have been pro- 
posed [H-HI whose large-distance dynamics is expected to 
be very different from the one of QCD. In particular it 
was shown that models falling in these frameworks could 
satisfy the experimental constraints (for recent reviews 
of techicolor models see [ltjjllj]). 
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Good candidate models in these frameworks are those 
which lie close to the lower boundary of the so-called 
conformal window, where the presence of an (approxi- 
mate) IR fixed point is believed to significantly change 
the non-perturbative dynamics of the theory. 

The use of matter fields in higher dimensional repre- 
sentations has been recently advocated [3, EH as an ef- 
fective and economic way of satisfying all known experi- 
mental constraints. 

In this work we will focus on one of these candidate the- 
ories, the so-called Minimal Walking Technicolor theory, 
based on the gauge group SU(2) with two Dirac fermions 
in the adjoint representation. 

All viable candidate models share the common prop- 
erty of being strongly interacting at the electroweak scale 
and as such are not fully under control by analytic meth- 
ods only. Although the analytical approaches are indis- 
pensable to show which models are the most promising, 
they all depend on uncontrolled approximations or con- 
jectures based on educated guesses. 

In this work we study MWT using first-principle nu- 
merical simulations, which allow to investigate the full 
non-perturbative dynamics of the theory. We use the 
same techniques matured during the last decades for Lat- 
tice QCD and which are now source of reliable and valu- 
able informations for the phenomenology of the strong 
interactions at high-energy experiments. Within the Lat- 
tice Gauge Theory framework quantitative predictions 
can be obtained, which demonstrate if a candidate model 
is indeed viable or not. 

In the last two years renewed interest among the lattice 
community has led to an increasing number of studies by 
several different groups [l6l - [55l |. 

In this work we present a detailed study of the non- 
singlet mesonic sector of the spectrum of the gauge theory 
SU(2) with two Dirac adjoint fermions. In a companion 
paper [56( we will present our result for the glueball mass 
spectrum and string tension and compared them to the 
one obtained in this paper for the mesonic spectrum as 
first suggested in [42[ . As the simplest of such interesting 
models, it is particularly amenable to numerical investi- 
gations. 
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Given the present analytical uncertainties, it is not 
clear if this theory lies within the conformal window or 
not. To understand if this model lies within the confor- 
mal window, in this paper we will compare our data to 
the signatures of spontaneous chiral symmetry breaking 
on the one hand and to the expected scaling behavior in 
proximity of an IR fixed point on the other. By study- 
ing the dependence of the low-lying meson masses on the 
current quark mass we will provide evidence for the ex- 
istence of an IR fixed point. 

This paper is organized as follows. In Sect. [TT] we re- 
mind to the reader the physical implications of the exis- 
tence of an IR fixed point in the theory and its observable 
consequences as derived from a Renormalization Group 
(RG) analysis of the fixed point. In Scct. lIIII we introduce 
the Lattice Gauge Theory formalism used this work. In 
Sect. IIVI we present our numerical results and compare 
them to the theoretical expectations. We finally conclude 
in Sect. El 



II. NON-QCD BEHAVIOR 

Strongly interacting theories that have a dynamics dif- 
ferent from QCD are needed in order to be able to build 
successful phenomenological models of DEWSB. One ex- 
ample of such theories is provided by gauge theories in 
the so-called conformal window, which are characterized 
by the existence of an infrared fixed point in their renor- 
malization group flow. 

A four-dimensional gauge theory minimally coupled to 
fcrmions in some representation of the gauge group has 
a perturbative UV fixed point provided the number of 
fermion species is not too large; at the UV fixed point 
the gauge coupling g and the fermion mass m are rele- 
vant couplings. The gauge coupling g is dimensionless 
in four dimensions. It is convenient to use dimension- 
less couplings for discussing RG flows. Hence we shall 
consider the dimensionless quantity rh = am = m//i 
when studying the RG transformations of the fermion 
mass. Renormalizcd trajectories, i.e. lines of constant 
physics, are one-dimensional curves originating from the 
UV fixed point. Points on a given renormalized trajec- 
tory correspond to theories that have the same long- 
distance physics, but different values of the UV cutoff. 
Each line corresponds to a theory with a given physical 
fermion mass. 

One of the lines of constant physics corresponds to the 
massless renormalizcd trajectory. If the theory possesses 
an IRFP, the latter has to lie on the massless trajectory, 
otherwise the finite fermion mass would drive the theory 
away from the IRFP at large distances. Assuming the 
existence of such a fixed point, we can linearize the RG 
equations in the vicinity of the fixed point, and identify 
relevant and irrelevant directions. In particular the mass 
rh will be a relevant operator. 

The running of the couplings is described by the RG 



equations: 



/i-^-m = -j(g)m, 
djji 



(1) 

(2) 



where g and m are respectively the running coupling and 
the running mass, which depend on the energy scale /i. 
Note that chiral symmetry guarantees that the RHS of 
Eq. @ is proportional to the mass itself. The function 7 
is the anomalous dimension of the scalar density operator 
-ip{x)ip{x). Eq. © implies: 



d/i \ /i 



[7+1] 



(3) 



Similar equations describe the evolution of all the other 
couplings that are compatible with the symmetries of the 
system under study. We shall denote the generic, di- 
mensionless coupling gi] their evolution is dictated by a 
corresponding /3 function: 



(4) 



The IRFP is defined by an isolated zero of the (3 func- 
tions. Theories in the conformal window become scale- 
invariant at large distances, and therefore cannot develop 
condensates. In particular chiral symmetry cannot be 
spontaneously broken, and there are no single-particle 
states; the dynamics is entirely expressed by the expo- 
nents that characterize the power-law behavior of field 
correlators at large distances. 

In the vicinity of a fixed point g* the RG equations 
can be linearized; the evolution is characterized by the 
matrix: 



Rij — 



9gj 



(5) 



The evolution of the dimensionless eigenvectors of the 
matrix R, Ui, is given by simple power laws: 



Ui(n) tx /i Vi 



(6) 



where yi are the eigenvalues of R^. The yi are the crit- 
ical exponents that are commonly used in the theory of 
critical phenomena. It is clear from Eq. ^ that yi > 
characterizes the relevant directions at the IRFP. 

The fermion mass is a relevant operator at the fixed 
point, and we can readily deduce from Eq. (J3j) : 



y m — 7* + l 



(7) 



where 7* is the value of the anomalous dimension at the 
fixed point. 

The scaling (or conformal) dimension of the scalar den- 
sity A m is related to the critical exponent y m by: 



Vn 



D 



(8) 
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where D is the dimension of space-time. The scaling 
dimension for a scalar operator is bound to be greater 
than one by unitarity, and it is equal to three for the 
scalar density in the free theory. This corresponds to the 
usual range < 7* < 2. 

Scaling of the free energy density. Let us now 
consider a RG transformation such that the lengths are 
rescaled by a factor b: 

a' = ba, [i' = (i/b. (9) 

Following the discussion above, the transformation of the 
singular part of the free energy density under such trans- 
formation can be written as: 

f s ( Ui , am, a/L) = b- D f s (b y > Ul , b y ™am, ba/L) . (10) 

We have denoted by Ui the irrelevant operators at the 



IRFP, and therefore yi < 0, while y m > as expected 
from the discussion above on the role of the fermion mass. 
We have included the dependence on the size of the box 
L, since this will provide finite-size scaling laws. The 
inverse of the size 1/L is treated as relevant coupling 
with unit eigenvalue; the underlying hypothesis here is 
that the finite value of L does not affect the RGE for 
the other couplings, i.e. that it is larger than the inverse 
mass of the states in the theory, Lm <C 1. 

Iterating the RG transformation n times yields: 



f s (ui, am, a/L) = f s {b ny *u l: b ny ™am, b n a/L) . 

(11) 

Choosing n such that b nVm am = amo, where mo is some 
reference mass scale, Eq. (jlip can be rewritten as: 



D / Vj 1_ \ D I \vj \ 

f m\ Vm ( ( m \ Vm ( m \ <>m \ f m\ ym I , / m \ vm \ 

f s{ u l ,am,a/L)={-) f s ^- j = (- j (- j ^.(12) 



Eq. (fT2l) describes the scaling with the fermion mass, and 
the functional form of finite-size effects. Expanding $ in 

/ \\vi\/y m 

powers of = ( J Ui, assuming as usual that f s 

is analytic as a function of the irrelevant couplings, yields 
the corrections to the scaling. Note that these corrections 
vanish when x = 0. 

Ignoring the corrections to scaling we obtain an expres- 
sion for the free energy that is useful to derive finite-size 
scaling properties: 



f s (am,a/L) = — 



D/y„ 



<S>{mL ym ). (13) 



Scaling of correlators. 

lators: 



Two-point function corre- 



f H {x\am) 



jD-l 



x(H(x)H(0) 1i ) 



(14) 



satisfy similar RG equations: 



d d 

a- 7m- 27^ 

oa dm 



/ fl (i;am)=0, (15) 



Iterating this relation n times, and performing manipu- 
lations analogous to the ones above, yields: 



fnix; am) 



\m J 



2"1H /Vr, 



* (x/|m/m r 1/am ) . (17) 



If the correlator decays exponentially at large distances 
with the mass of the lightest state that overlap with the 
fields, then Eq. (|17l) determines the scaling of this mass: 



M - m 1/y " 



(18) 



Using the relation between y m and 7* discussed above, 
the scaling for the masses of the physical states is 



M 



,1/(1+7.) 



(19) 



Note that according to this analysis all masses scale with 
the same exponent. Different states are selected by dif- 
ferent fields appearing in the correlator which implies 
that the anomalous dimension in the prefactor, yn, does 
change. However the scaling of the mass is entirely due 
to the RG behavior of the argument of the function ^, 
which does not depend on the state under scrutiny. 



where we have neglected the dependence on the irrele- 
vant couplings, and 7# is the anomalous dimension of 
the field H . The solution to the above equation obeys 
the following scaling law: 



f H (x;am) = b~ 2 ^ H G(x;b y ™ am) 



(16) 



III. LATTICE SIMULATIONS 

Non-perturbative numerical simulations are performed 
after introducing an effective UV and IR cutoff in the 
form of a space-time lattice of finite extent. 
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The Euclidean path integral is thus reduced to an or- 
dinary integral over a large number of degrees of free- 
dom. The choice of a discretized action on the lattice 
is not unique; different choices will result in different 
lattice artefacts. Unimproved Wilson fermions are used 
throughout this work. The lattice action with matter 
fields in a representation R is given by 1 

"/ 

S(U,^^) = S g (U) + Y,^D m (U R )A , (20) 
»=i 

with explicit expressions for the gauge action S g and the 
massive lattice Dirac operator D m given below. For this 
choice of discretization, the action depends only on two 
bare parameters: the bare inverse coupling /3 and the 
bare dimensionless quark mass amo, a being the lat- 
tice spacing. While the link variables appearing in the 
gauge action are in the fundamental representation 2 of 
the gauge group, the links in the lattice Dirac operator 
are in the same representation R as the fermion fields. 
The partition function, after integrating out the matter 
fields, takes the form: 

Z = J cxp[-S g (U)][dct D m {U R )] n fdU . (21) 

For the Wilson action used in this work the gauge ac- 
tion is given by 

S i( U ) = Jf E Re tr V ^ x ) > ( 22 ) 

where V^vix) is the elementary lxl plaquette in the fj,— v 
plane at lattice site x. In the fermion action the Wilson- 
Dirac operator is given by 

D m (U R ) = am + ~ £ r> (V m + V£) - aV^V M ] , 

(23) 

where is the discretized forward covariant derivative 
depending on the link U R and V* its adjoint operator: 

(V^)(x) = U R (x, »)tP(x + - i>{x) . (24) 

In a numerical lattice simulation the computation of 
the discretized path integral is performed by Monte Carlo 
integration using importance sampling: an ensemble of 
gauge configurations is generated with probability pro- 
portional to exp[-0S g (U)] [dct D m (U R )] n f. The expec- 
tation value of any observable can then be computed as 
a stochastic average over this ensemble of configurations. 



1 We omit for the sake of simplicity all the position, color and spin 
indexes. 

2 The link variables can also be taken in a different representation 
of the gauge group if one chooses to. 



A. Sources of systematic errors 

In order to obtain continuum values for the observables 
of a theory from numerical simulations of its lattice dis- 
cretised version, an appropriate limiting procedure must 
be performed. It is thus important to understand when 
and to what extent the outcome of lattice simulations 
are a faithful depiction of the continuum physics. This is 
especially important when one tries to understand a new 
theory as the MWT in the present work, since we lack the 
insight and the experimental input we have for instance 
in the more familiar case of QCD. In fact, in order to be 
the description of a new force of Nature, this theory has 
to be rather different from QCD and we need to ensure 
that we are observing genuine features of the continuum 
theory, and not artefacts of our lattice formulation. 

We will now list the most important sources of sys- 
tematic errors which are present in lattice simulations 
and what are the appropriate limits to take in order to 
recover the continuum physics. In addition to these, sta- 
tistical errors are also always present, but those can be 
reduced arbitrarily by producing a big enough ensemble 
of configurations. 

Finite-size, finite-temperature effects. These 
are due to the presence of an IR cutoff in the form of a 
finite extent of the 4-dimensional lattice both in the spa- 
cial and temporal directions. The standard lattice geom- 
etry used in numerical simulations is TxL 3 , i.e. the three 
spacial directions have equal length. The correct vacuum 
expectation values of the continuum theory are recovered 
in the limit in which T, L — > oo. On a 4-dimensional torus 
these expectation values can have large corrections, even 
if asymptotically the infinite volume limit is reached at 
an exponentially fast rate. As the system is tuned closer 
to a critical point, the magnitude of finite-size effects 
and the autocorrelation of lattice observables increase. 
This is what happens for example when particles with 
a Compton wavelength comparable with the lattice size 
are present. Moreover, if the size of the lattice is not suf- 
ficiently large, the system may enter a phase that bears 
little resemblance to the large-volume theory we are in- 
terested in. Some examples, based on analytical finite 
volume results (57rl59| . have been recently discussed in 
Ref. El- 

Explicit breaking of chiral symmetry. It is dif- 
ficult and numerically very expensive to preserve chiral 
symmetry once the theory is discretized on a lattice. For 
the particular choice of the lattice action used in the 
work, i.e. Wilson-Dirac fermions, chiral symmetry is 
explicitly broken at the Lagrangian level. In order to 
recover it in the continuum limit the tuning of one pa- 
rameter (the bare quark mass) is necessary. Even with 
lattice chiral fermions, the low modes of the lattice Dirac 
operator, which appear as the chiral limit is approached, 
make it practically impossible to run simulations at very 
light masses with the current algorithms. When using 
Wilson fermions on a given finite lattice, there is a lower 
limit to the masses that are numerically accessible. This 
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limit depends on the volume of the system [6CJ, |6l| . 

Discretization artefacts. The space-time lattice 
introduces a UV cutoff, i.e. the lattice spacing a. The 
interesting continuum physics is recovered in the limit 
where the lattice spacing goes to zero. However since all 
the quantities in a simulation are dimensionless, the lat- 
tice spacing a is not a parameter that one can directly 
change. Instead the value of the lattice spacing is in- 
ferred from the measure of a physical quantity on the 
lattice, which is thus used to fix the scale a. The con- 
tinuum limit is recovered when the system is tuned to 
a UV critical point where the lattice spacing vanishes. 
This can be done in the space of bare parameters by 
increasing the lattice inverse coupling f3 to infinity, i.e. 
sending the coupling of the theory g to zero, as dictated 
by asymptotic freedom. At that point the system under- 
goes a continuous phase transition and, by universality, 
the microscopic details of the system become immaterial. 
To put it in a different way, the physics at the scale of 
the lattice spacing must decouple from the long range 
(continuum) physics, which is what we are interested in 
studying on the lattice. 

To give reliable predictions, lattice simulations must 
be performed in a regime in which a precise hierarchy 
of scales is realized. To illustrate this point, let us con- 
sider first the case of QCD. Numerical simulations, as ex- 
plained above, are always performed with a finite quark 
mass, which explicitly induces a mass gap in the theory. 
To avoid finite-size effects in the computation of the low- 
lying spectrum, the lattice size L must be much bigger 
than the inverse mass of the lightest hadron Mps we aim 
to measure. If we are interested in the chiral regime, 
as it is usually the case, this light hadron mass must be 
much smaller than the characteristic hadronic scale at 
which the theory becomes strongly interacting. In QCD 
we can take for example the Sommer scale ro as a con- 
venient quantity to measure on the lattice. Finally to 
avoid discretization errors, the lattice spacing must be 
much smaller than the reference scale for strong inter- 
actions ro, so that the physics at the scale of the UV 
cutoff is weakly interacting. In summary we must have 
the following hierarchy of scales: 

L- 1 < M PS < r^ 1 < a' 1 , (25) 

for the computations to reproduce reliably the features 
of the continuum chiral regime of QCD. 

Let us now consider the case in which the underlying 
continuum theory has an IR fixed point in the massless 
limit. The presence of a mass term in numerical simula- 
tions explicitly breaks conformal invariance and as in the 
QCD-like case a mass gap and a particle spectrum for 
the theory are generated. In order to reliably estimate a 
hadron mass, we still need the lattice to be big enough to 
fit the state we are interested to measure, i.e. as above 
<C Mps, for example. However in contrast to the 
QCD-like case, all the masses of the hadrons vanish as 
we send the explicit symmetry breaking term to zero. 
One can still define an IR scale like ro but, as the theory 



is no longer confining, ro is no longer related to the par- 
ticle masses and it is therefore not a useful quantity to 
compare to. One can be tempted however to introduce a 
modified version of the IR scale A which marks the onset 
of the IR scaling region. With this definition the interest- 
ing mass region to explore is Mps <C A, in analogy with 
the previous case. Finally also in this case discretiza- 
tion artefacts should be suppressed using fine lattices. In 
summary we have a similar hierarchy of scales as before 

L- 1 < M PS «A« a -1 , (26) 

with a new IR scale A whose definition is related to the 
new features of the theory. Unfortunately a convenient 
definition of A readily measurable on the lattice, like the 
Sommer scale for QCD, is not at hand. As defined above 
A can only be inferred a posteriori by observing the scal- 
ing behavior of some physical observable. 

B. Simulation Code 

The results presented in this work are obtained using 
our own simulation code, written from scratch for the 
specific purpose of studying gauge theories with fermions 
in arbitrary representations. The code was developed 
during the last few years and has been presented and 
tested in detail in Ref. [2lj . This code was designed to 
be flexible and easily accommodate for fermions in arbi- 
trary representations of the gauge group, without com- 
promizing the performance and ease of use of the code 
itself. 

Our simulation code, named HiRep, is suitable for the 
study of gauge theories with 

• gauge group SU(N) for any N. The code has 
already been used for the study of the large-N 
mesonic spectrum [62| with TV up to 6. 

• generic fermion representations. At present the 
code implements the fundamental (fund), adjoint 
(AD J), 2-index symmetric (S) and antisymmetric 
(AS). All of these different representation have al- 
ready been successfully tested and used [63j . It is 
easy to extend the code to other representations, 
like 3-index symmetric for example, and even to 
have fermions in two or more different represen- 
tations at the same time. In particular no mod- 
ifications for the computation of the HMC force 
are needed, which is typically the most complicated 
part of the code. 

• any number of flavors. We use Wilson fermions 
with the HMC 64]/RHMC[6l[6l] algorithm, which 
is an exact algorithm for any number of flavors. 

We have also implemented a significant number of observ- 
ables such as the the measure of the mesonic spectrum 
(presented in this work), Schrodinger functional observ- 
ables (used in Ref. 48]), gluonic observables like Wilson 
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and Polyakov loops, glueball masses presented in a com- 
panion paper (5(| . 

Since this was a new code, and lattice simulations in 
the past were mainly devoted only to QCD, we made a 
large effort to validate the code and to study its behavior 
in the parameter region of interest for the physics. As 
part of this effort, in particular we made a number of 
cross-checks: 

• for SU(3), rif = 2 with different, established codes 
(we used M. Luscher's DD-HMC algorithm 

• consistency among different representations (e.g. 
SU(3) AS vs fund, SU(2) SYM vs ADJ); 

• observables with different codes (e.g. meson spec- 
trum for SU{2) ADJ in Chroma); 

• large quark mass limit compared to pure gauge 
(quenched) spectrum; 

• correctness of integrator, reversibility, acceptance 
probability; 

• independence from integrator step size. 

To control the stability of the simulations, which 
could incur in well-known problems close to the chiral 
limit [6(| , we monitored the lowest eigenvalue of the (pre- 
conditioned) Wilson-Dirac matrix. The average value of 
the lowest eigenvalue together with the standard devia- 
tion of their distribution can be found in the tables in 
Appendix [Bl We found no instabilities in our runs. 

C. Meson masses 

Let r and V be two generic matrices in the Clifford 
algebra, we define the two-point correlator at zero mo- 
mentum as follows: 

frr>(t) =^(Wi(x,i)rf 2 (x,i))tf(0)^ 2 (0)), (27) 

X 

where "01 and ip2 represent two different flavors of de- 
generate fermion fields, so that we only consider flavor 
non-singlet bilinears. Denoting the space-time position 
(x, t) by x and performing the Wick contractions yields: 

fvr- (*) = J2 - tr [7 r t 7o 5(.T, 0)r'S(0, x)] , (28) 

X 

where S denotes the quark propagator, i.e. the inverse 
of the hermitean Wilson-Dirac matrix 75-D. In practice 
not all matrix elements of S are computed, but only some 
single rows (point-to-all propagator) by solving the linear 
system: 

D(x, y)ABVB'°(y) = S A,A S x,o , (29) 

where capital Latin letters like A = {a, a} are collective 
indices for color and spin, and A, x = is the position 



of the source for the inverter. The inversion is performed 
using a QMR recursive algorithm with even-odd precon- 
ditioning of the Dirac operator, which is stopped when 
the residue is less than 10 -8 . For some of our lattices we 
used the noise-reduction technique described in Ref. [68[ 
to take a stocastic average over the volume of the point 
source. 

Following Ref. (6f[, masses and decay constants for 
the pseudoscalar meson are extracted from the asymp- 
totic behavior of the correlators /pp and /ap at large 
Euclidean time. The pseudoscalar mass and the vacuum- 
to-meson matrix element are obtained from the correla- 
tor of two pseudoscalar densities: 

/pp(t) = exp [-Mp S f] + • • • ■ (30) 

The meson mass is obtained by fitting the effective mass 
to a constant, while the coupling Gps is extracted from 
the amplitude of the two-point function /pp. The defi- 
nition of the effective mass used in this work is given in 
Appendix |A] 

As in Ref. 61] the ratio 

m cS (t) = \ \(d + d*) / AP (t)] // PP (t) (31) 

yields the PC AC mass m with corrections of 0(a) for 
the unimproved theory. Note that the decay constant 
is not computed directly; it is obtained from the values 
computed above as: 

^ p s = T72~ G ps- (32) 

The decay constant extracted from bare lattice correla- 
tors is related to its continuum counterpart by the renor- 
malization constant Za, which has been computed in per- 
turbation theory in Ref. [1 8j . 

Finally the mass of the vector state is extracted from 
the /w correlator, again using a fit to the effective mass 
plateaux. 

On the smallest lattices that we have used in this study, 
it is difficult to isolate clearly the contribution from the 
lowest state, which dominates the large-time behavior of 
two-point correlators; this yelds large systematic errors. 
We have however explicitly measured the meson masses 
at the same value of the bare parameters - corresponding 
to the same value of the PC AC mass - on increasingly 
larger lattices to ensure that the residual finite-volume 
corections on these quantities are small. 

For the pseudoscalar and vector channels we have also 
verified that different choices of interpolating operators 
give results that agree with each other. In these channels 
the final results presented below are obtained taking the 
average over the different choices of interpolating opera- 
tors. 
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D. Simulation parameters 

All the simulations discussed in this work are per- 
formed at a fixed lattice spacing, corresponding to a bare 
coupling (3 = 2.25. This value of the coupling was chosen 
based on previous studies of the same theory 0, [22|, |34| 
to avoid a bulk phase transition present at about (3 = 2.0. 
The extrapolation towards the continuum limit requires a 
new series of runs at different values of the bare coupling, 
and is left to future investigations. 

We use four different lattices: 16 x 8 3 , 24 x 12 3 , 32 x 16 3 
and 64 x 24 3 . For each of these four lattices a num- 
ber of ensembles corresponding to different quark masses 
were generated, focusing in particular on the range cor- 
responding to pseudoscalar masses between 0.6a _1 and 
0.2a -1 . As explained above, when using Wilson fermions, 
the chiral and infinite volume limits are intertwined. In 
practice for the simulation to be stable, one cannot arbi- 
trarily decrease the quark mass without also increasing 
the volume. This is also necessary to keep under control 
finite- volume effects, thus remaining in the large volume 
limit. We explicitly control the size of these systematic 
errors performing simulations with different volumes. 

For each lattice and quark mass we accumulated a sta- 
tistical ensemble of about 5000 thermalized configura- 
tions, except at the largest volume for which we present 
only preliminary data based on approximately 500 con- 
figurations for each quark mass. The gauge configura- 
tions were generated using trajectories in the molecular 
dynamics integration of length 1 for the two smallest vol- 
umes and 1.5 for the two largest lattices, with integration 
parameters leading to an acceptance rate of about 85% 
in all cases, and to an integrated autocorrelation times 
for the lowest eigenvalue of the Wilson-Dirac operator of 
order 15 or less. 

Details of simulation parameters and results are re- 
ported in the tables of Appendix [B] 

IV. SU(2) WITH 2 ADJOINT FERMIONS 



chiral limit the pseudo-scalar particles (pions) become 
massless, while the other states in the spectrum remain 
massive. In the small mass regime the theory can be 
effectively described by a chiral lagriangian and the fa- 
miliar results of QCD can be recovered. In particular the 
theory has a non-zero pseudoscalar decay constant Fp§ 
and chiral condensate (ipi/)). Using the PCAC mass from 
the axial Ward identity, to, to parametrize the explicit 
breaking of chiral symmetry, we expect the usual scaling 
Mp S oc to, as to — y 0. The Gcll-Mann-Oakes-Renner 
relation is satisfied and can be used to extract the chiral 
condensate: (Mps-Fps) 2 /m — >• —(tpip). 

On the other hand, if the theory has an IR fixed point, 
the arguments of Sect. [TT| apply. In the scaling region, i.e. 
in proximity of the IR fixed point, dimensionful phys- 
ical quantities are expected to scale with a power law 
behavior. The universal exponents appearing in these 
scaling laws are related to the anomalous dimensions of 
the scaling fields as shown in Sect. Ql] In particular the 
scaling of the hadron masses is governed by the anoma- 
lous dimension of the mass 7* at the fixed point. If we 
parametrize the explicit breaking of chiral and conformal 
symmetry by to, in the massless limit m — > we expect 
that all hadron masses vanish proportionally to the same 
power of to: Mh a d ~ m 1 ^ 1 " 1 " 7 *); in particular the ratio 
of the vector to pseudoscalar meson masses remains fi- 
nite: My /Mps — > const < 00. Also the pseudoscalar 
decay constant Fpg and the chiral condensate {ipip) are 
expected to vanish in the chiral limit. 

The behavior just described assumes that the system is 
in an infinite volume near the continuum limit. However 
this is not always the case in a lattice simulation. In 
particular the finite size of the system can be seen as a 
relevant coupling which will drive the system away from 
criticality under the renormalization group (RG) flow. 
As to is decreased towards the chiral limit, great care 
must be taken to control finite size effects as explained 
in Section Unl We will show below that these effects are 
quite big in the range of masses and lattice volumes which 
are currently used in numerical simulations. 



Before looking at the actual numerical results, let us 
discuss the signatures of an IR fixed point, in order to 
focus on the important aspects of our numerical evidence. 
In this work, we search for indications of IR conformal 
behavior in the spectrum. This is not the only possible 
way, as one can, for example, study the non-perturbative 
running of a coupling defined in some particular scheme. 
In fact most of the claims of the existence of IR fixed 
points so far have been made by looking at the evolution 
of the c oup ling in the Schrodinger functional scheme [l?], 

MM 

UK]. However at present these studies still lack 
a reliable continuum extrapolation and the claims of the 
existence of IR fixed points should therefore be confirmed 
by more solid numerical investigations. 

If there is no IR fixed point, the theory is expected 
to be confining and chiral symmetry to be spontaneously 
broken. We will refer to this case as QCD-like: in the 



A. Spectrum 

For each ensemble of configurations, we measured the 
quark mass from the axial Ward identity (PCAC mass) 
to, the pseudoscalar meson mass -Mps, the vector mass 
My and the pseudoscalar decay constant Fps- 

We locate the chiral limit at the critical bare mass 
where the PCAC mass vanishes. This does not corre- 
spond to zero bare mass because the explicit breaking 
of chiral symmetry with Wilson fermions induces an ad- 
ditive renormalization of the quark mass. We show in 
Fig. [1] the extrapolation of to for different lattice sizes. 
Using a linear extrapolation of the four lightest measured 
points, the chiral limit can be located at the critical bare 
mass am c — —1.202(1). As expected from the fact that 
to is an UV quantity, no significant finite-size effects are 
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FIG. 1: Extrapolation of the quark mass from the axial Ward 
identity to locate the chiral limit. As expected no significant 
finite size effects are present. 



visible and the measured values for this quantity agree 
within errors on all four lattices. 

Our results for the mass of the pseudoscalar meson 
are presented in Fig. [2] The interesting region of small 
quark masses is shown in the right panel. Given the 
level of accuracy of the present measure, the finite volume 
systematics on Mpg are clearly visible and, as the PCAC 
mass is decreased, they become more and more relevant, 
as discussed in Sect. IIII1 To quantify this systematic 
effect and to keep it under control, we use larger lattices 
as the chiral limit is approached. 

These large finite size effects make it harder to draw 
definitive conclusions about the functional behavior of 
the pseudoscalar mass in the chiral limit. For a QCD- 
like theory the ratio aMp S /m, shown in Fig- El should be 
a (non-zero) constant in the chiral limit. On the other 
hand if the theory has an IR fixed point the ratio should 
vanish in the chiral limit if 7* < I or diverge if 7* > f . 
Our data clearly favor the IR conformal scenario with 
7* < 1. An accurate determination of the anomalous 
dimension is more difficult, but from the almost linear 
behavior of Mp$ as a function of m a small value of 7* 
seems to be preferred. 

The IR conformal scenario is also favored when one 
looks at the ratio of the vector to the pseudoscalar mass 
which is shown in Fig. [4] This quantity is bounded to be 
greater than 1 [6^| and in the heavy quark limit will ap- 
proach unity. At large m finite volume effects are small, 
the ratio is bigger than 1 and decreasing as m increases, 
as expected in the heavy quark approximation. What 
is remarkable in our data is the fact that in the whole 
mass range we were able to explore, and in which the 
pseudoscalar mass changes roughly by a factor of 7, the 
vector meson never becomes more than 5% heavier than 
the pseudoscalar, so that the ratio remains approximately 
constant in the chiral limit. This is the expected behav- 
ior in an IR conformal theory, since in this case all the 



hadronic masses scale with the same critical exponent. 

Another physically interesting quantity to consider is 
the pseudoscalar decay constant Fps, shown in Fig. [5] 
Among the ones presented in this paper, this is the quan- 
tity which shows the largest sensitivity to finite- volume 
effects. By looking at the behavior of Fps at different 
volumes, an envelope of the curves as a function of the 
PCAC mass is clearly visible, which should be used for 
the chiral extrapolation. For a QCD-like theory the re- 
sult in the chiral limit is a non-zero value. The direct 
extrapolation however is difficult to carry out with rea- 
sonable accuracy; for this reason, we prefer to exploit the 
finite-size effects themselves to obtain a more insightful 
statement. As discussed in Sect.[II]near an IR fixed point 
one can consider the finite size L of the system as a rele- 
vant parameter in the RG flux and thus obtain universal 
scaling laws for physical observables, Eq. (IT3l) . This fi- 
nite size-scaling law can be conveniently rewritten, for 
example for the pseudoscalar decay constant, as: 

LFp S = T(W/( 1+7 ')). (33) 

Scaling is observed if the different curves corresponding 
to keeping the volumes fixed and varying the quark mass, 
collapse on top of each other. As a byproduct of the 
procedure an estimate of the critical exponent is also ob- 
tained. To illustrate the procedure, we plot LFps as a 
function of x = Lm 1 ^ 1+1 *^ in Fig. [5] for various values 
of 7*. Good scaling is observed for 0.05 < 7* < 0.20, 
while larger values of 7* (and in particular 7* = I) seem 
to be excluded. The observed scaling is again in agree- 
ment with the existence of an IR fixed point with a small 
7* in the MWT theory. The range of values of 7* for 
which a good quality of the scaling is obtained is com- 
patible with independent estimates performed with the 
Schrodinger functional [48| . 

A phenomenologically relevant quantity to look at is 
the ratio My/Fpg- This ratio is shown in Fig. [7l The 
large dependence of Fps on the finite-size of the lattice 
is reflected in the large finite-size effects for the ratio. A 
tentative large- volume limit curve can be obtained by dis- 
carding the results at the lightest masses on the smaller 
volumes. As the fermion mass drops below 0.1a -1 , the 
ratio starts decreasing unless the volume is made larger. 
Taking the envelope of the curves for different values of 
the volume, one can expect a value of about 5 — 6 in the 
chiral limit. 

The chiral condensate would also be a prime candidate 
to study chiral symmetry breaking. However due to the 
use of Wilson fermions, the direct measure of (^"0) is 
plagued with UV divergences which are notoriously dif- 
ficult to tame. Using the GMOR relation an estimate 
for the chiral condensate can be obtained 3 . The method 



3 We do not attempt here to compute the necessary multiplicative 
renormalization constant, since we are not interested to the ac- 
tual physical value. Perturbative results for the renormalization 
of fermions bilincars can be found in Ref. UM 
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FIG. 2: Pseudoscalar meson mass as a function of the PCAC mass. The interesting small mass region shaded in the left panel 
is enlarged on the right. Finite volume effects are evident and grow approaching the chiral limit. 
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FIG. 3: Ratio of the pseudoscalar mass squared to the PCAC 
mass. The extrapolation to the chiral limit suffers from large 
finite-volume effects. See the text for a discussion. 
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FIG. 4: Comparison between the vector and pseudoscalar me- 
son masses. At large PCAC mass, due to quenching the ratio 
is very near to one. Near the chiral limit large finite size 
effects show up. 



FIG. 5: Pseudoscalar decay constant near the chiral limit. 
Very large finite volume effects are present also in this case 
which cause the chiral extraplation to be have large uncer- 
tainties. 



has been applied with success in the case of QCD, see 
e.g. Ref. [7(|. We present our results for this quantity 
in Fig. HI Although there is a partial cancellation of the 
finite size effects coming from the pseudoscalar mass and 
the decay constant, the larger volume dependence of the 
latter dominates, yielding large systematic errors. As a 
consequence an extrapolation is unfortunately not possi- 
ble from our current set of data. We observe that finite 
volume effects tend to make the condensate smaller, how- 
ever the small numerical value of the bare condesate by 
itself is not meaningful: for example in a typical QCD 
simulation the value for this quantity is an order of mag- 
nitude smaller than the one presented here. 
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FIG. 7: Vector to pseudoscalar decay constant ratio. Large 
finite-size effects are present also in this case which make the 
extrapolation to the chiral limit difficult. The envelop of the 
curves in the plot suggests a limit value of about 5 — 6. 



FIG. 8: The GMOR relation can be used to extract informa- 
tion on the chiral condensate. The measure results however 
quite difficult in practice and we cannot distinguish any signal 
of spontaneous chiral symmetry breaking. 



V. CONCLUSIONS 

In this work we have presented a careful investigation 
of the mesonic spectrum of one of the candidate theo- 
ries for a realistic technicolor model, the so-called Min- 



imal Walking Technicolor, based on gauge group SU(2) 
with two Dirac adjoint fermions. Theoretical specula- 
tions about this theory indicate that it is very near to the 
lower boundary of the conformal window. In this work 
we used numerical lattice simulations to look at mesonic 
spectrum and we found some evidence that the theory 
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lies in fact inside the conformal window and possesses an 
IR conformal fixed point. 

Such numerical simulations are an extremely powerful 
tool to explore the non-perturbative dynamics of gauge 
theories which is otherwise inaccessible to theoretical 
speculations, but great care must be taken to control sys- 
tematic errors. In order to tame finite size corrections, 
which make the extrapolation to the chiral limit difficult, 
in this work we aimed for the first time at reaching the 
chiral limit in a controlled way: we used a series of four 
different lattice sizes up to a large 64 x 24 3 . 

Evidence for the existence of an IR fixed point was 
found in the behavior of the different mesonic observ- 
ables analyzed, namely the pseudoscalar and vector me- 
son mass and the pseudoscalar decay constant, which 
show significant deviations from the expectations of a 
more familiar QCD-like scenario, where spontaneous chi- 
ral symmetry breaking occurs. We showed that our data 
are compatible with the existence of an IR fixed point by 
using the predicted scaling laws that need to hold in this 
case. 

Although the present data show clear signs of confor- 
mality in the infrared, our study still has several limita- 
tions which should be addressed in the future to put our 
results on a more solid ground. Smaller quark masses 
and consequently larger lattice volumes would increase 
the reliability of the scaling analysis we performed in this 
paper. However the major source of uncertainty is the 
fact that all numerical simulations used in this work were 
performed at a single value of the lattice spacing, and no 
test to assure the validity of our findings in the contin- 
uum limit has been done so far. 

Finally in this paper we focused our attention only on 
the mesonic spectrum, while substantially more informa- 
tion can be gained by combining it with observables from 
other sectors of the theory, as we proposed in Ref. (42J. 
The detailed study of gluonic observables, and their com- 
parison to the mesonic ones is the subject of a companion 
paper [56| , which provides further evidence for the exis- 
tence of an IR fixed point. 
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Appendix A: Effective mass definition 

For the definition of the effective mass used in this work 
we follow Ref. 71]. A mesonic correlator on the lattice 
has the form: 



A I 



C(t) = a m cosh [E m T\ 



(Al) 



with t = t - T/2 = 0, 1, 2, ... , T/2, where we consider 
only M excited states. Now since: 

(cosh [E m ]) n = — ( k ) c °sh [E m (2k - n)} , (A2) 
fe=0 ^ ' 

taking similar linear combinations of the C(r) we have: 



1 n ( \ M 

J °( 2k - n ) = E a ™( cosh m 



2™ ^ \k 

k=0 



Introducing the variables 

x m = cosh [E. { 

1 - 
= 2™ E V /■ 

k=0 



mj ; 

n 



C(2k - n) 



we can rewrite Eq. IA3I in matrix form as 



(A3) 

(A4) 
(A5) 




'1111 

Xi X 2 X 3 Xi 
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V 



(A6) 



We need to solve Eq. lA6l where y n is known and both a m 
and unknown. It is always possible to find the 

unique solution to Eq. IA6I considering 2M consecutive 
points in the transformed correlator y n . In general the 
x n are given by the roots of the M-degree polynomial: 



det 



( Vo 2/i • • • Vm-i 1 \ 

Vl V2 ■■■ VM X 

2/2 2/3 ■■■ 2/Af+l x 2 

\UM 2/M+l ' ' ' V2M-1 X M J 



(A7) 
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and the a m are then given by the solution of the linear 
system Eq.(4) obtained with the known x n . 

In this work we do not consider excited states and we 
only need the solution of Eg. 1X71 for M = 1 which is given 
by x = yi/y . 



Appendix B: Tables 
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lattice 


V 


— OWio 




(P) 


T 


A 




AO 


16 x 8 3 


0.95 


7601 


0.63577(16) 


5.45(72) 


3.582(13) 


8.6(1.4) 


Al 


16 x 8 3 


0.975 


7701 


0.63843(15) 


5.43(71) 


2.982(12) 


6.65(96) 


A2 


16 x 8 3 


1 


7801 


0.64136(15) 


5.10(64) 


2.427(11) 


6.28(88) 


A3 


16 x 8 3 


1.025 


7801 


0.64463(15) 


4.29(50) 


1.894(10) 


6.07(84) 


A4 


16 x 8 3 


1.05 


7801 


0.64793(15) 


3.48(36) 


1.4596(79) 


4.39(52) 


A5 


16 x 8 3 


1.075 


6400 


0.65179(16) 


2.99(32) 


1.0692(74) 


4.27(55) 


A6 


16 x 8 3 


1.1 


6400 


0.65566(16) 


3.28(37) 


0.7564(60) 


3.77(45) 


A7 


16 x 8 3 


1.125 


7073 


0.66037(15) 


2.99(30) 


0.4854(43) 


3.03(31) 


A8 


16 x 8 3 


1.15 


6400 


0.66550(16) 


3.31(37) 


0.2779(31) 


2.80(29) 


A9 


16 x 8 3 


1.175 


6400 


0.67177(17) 


3.24(36) 


0.1351(18) 


2.80(29) 



TABLE I: Bare parameters and average plaquette for the 16 x 8 3 lattice. 



lattice 


-ar»o 


am 


aMps 


aM v 


aFps 


a 2 G PS 


AO 


0.95 


0.3899(40) 


1.4717(40) 


1.5203(52) 


0.3354(60) 


0.933(13) 


Al 


0.975 


0.3649(41) 


1.4093(43) 


1.4586(55) 


0.3240(61) 


0.883(13) 


A2 


1 


0.3365(42) 


1.3436(43) 


1.3936(58) 


0.3083(62) 


0.829(13) 


A3 


1.025 


0.3066(38) 


1.2630(46) 


1.3115(59) 


0.2908(58) 


0.759(12) 


A4 


1.05 


0.2749(39) 


1.1756(48) 


1.2233(64) 


0.2666(58) 


0.673(11) 


A5 


1.075 


0.2389(39) 


1.0623(56) 


1.1048(72) 


0.2392(59) 


0.567(11) 


A6 


1.1 


0.2031(39) 


0.9398(66) 


0.9784(86) 


0.2157(59) 


0.472(11) 


A7 


1.125 


0.1643(36) 


0.7817(76) 


0.811(10) 


0.1910(57) 


0.357(10) 


A8 


1.15 


0.1185(32) 


0.5740(89) 


0.587(11) 


0.1675(56) 


0.2347(82) 


A9 


1.175 


0.0650(24) 


0.330(11) 


0.3476(91) 


0.1611(50) 


0.1347(76) 



TABLE II: PCAC and meson masses from the 16 x 8 3 lattice. 



attice 


— amo 


am 


aMpg/m 


Mv/Fp S 


Mv/Mps 


a 3 (MpsFps) 2 /rn 


AO 


0.95 


0.3899(40) 


5.554(62) 


4.533(77) 


1.0330(12) 


0.625(19) 


Al 


0.975 


0.3649(41) 


5.443(67) 


4.502(79) 


1.0349(13) 


0.571(18) 


A2 


1 


0.3365(42) 


5.365(71) 


4.521(86) 


1.0372(15) 


0.510(17) 


A3 


1.025 


0.3066(38) 


5.203(70) 


4.511(84) 


1.0384(17) 


0.440(15) 


A4 


1.05 


0.2749(39) 


5.027(77) 


4.589(94) 


1.0405(22) 


0.357(13) 


A5 


1.075 


0.2389(39) 


4.724(83) 


4.62(10) 


1.0399(28) 


0.270(11) 


A6 


1.1 


0.2031(39) 


4.349(92) 


4.53(11) 


1.0410(37) 


0.2025(96) 


A7 


1.125 


0.1643(36) 


3.721(92) 


4.24(11) 


1.0375(53) 


0.1358(74) 


A8 


1.15 


0.1185(32) 


2.782(90) 


3.51(11) 


1.0242(88) 


0.0781(48) 


A9 


1.175 


0.0650(24) 


1.67(10) 


2.159(80) 


1.054(31) 


0.0434(30) 



TABLE III: Mass ratios from the 16 x 8 3 lattice. 
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lattice 


V 


-amo 




(P) 


r 


A 




BO 


24 x 12 3 


0.95 


10201 


0.635310(59) 


6.16(74) 


3.5058(50) 


3.08(26) 


Bl 


24 x 12 3 


1 


8652 


0.640998(64) 


4.92(58) 


2.4218(44) 


3.10(29) 


B2 


24 x 12 3 


1.05 


7819 


0.647633(70) 


6.79(99) 


1.4936(51) 


5.80(78) 


B3 


24 x 12 3 


1.075 


7186 


0.651630(68) 


4.61(58) 


1.0553(40) 


4.95(64) 


B4 


24 x 12 3 


1.1 


6393 


0.655827(76) 


4.09(51) 


0.7202(30) 


7.8(1.3) 


B5 


24 x 12 3 


1.125 


6200 


0.660588(75) 


3.97(50) 


0.4419(22) 


5.98(91) 


B6 


24 x 12 3 


1.15 


1599 


0.66588(15) 


3.71(90) 


0.2271(31) 


6.6(2.1) 


B7 


24 x 12 3 


1.175 


5582 


0.672074(79) 


4.22(58) 


0.08641(90) 


3.78(49) 


B8 


24 x 12 3 


1.18 


4081 


0.673474(92) 


4.01(63) 


0.06561(92) 


10(2.5) 


B9 


24 x 12 3 


1.185 


4201 


0.675094(93) 


3.42(49) 


0.05196(71) 


3.53(51) 


BIO 


24 x 12 3 


1.19 


3501 


0.67663(10) 


4.15(70) 


0.03985(61) 


5.2(1.0) 



TABLE IV: Bare parameters and average plaquette for the 24 x 12 3 lattice. 



lattice 


—amo 


am 


aMps 


aMv 


aFps 


a 2 G PS 


BO 


0.95 


0.3931(38) 


1.4746(23) 


1.5224(32) 


0.3343(62) 


0.925(12) 


Bl 


1 


0.3368(40) 


1.3495(26) 


1.4003(36) 


0.3020(63) 


0.819(12) 


B2 


1.05 


0.2765(40) 


1.1874(29) 


1.2383(40) 


0.2607(63) 


0.667(11) 


B3 


1.075 


0.2410(38) 


1.0809(30) 


1.1265(41) 


0.2320(58) 


0.5635(96) 


B4 


1.1 


0.2025(40) 


0.9614(35) 


1.0016(46) 


0.1991(53) 


0.4558(85) 


B5 


1.125 


0.1604(34) 


0.8020(41) 


0.8312(56) 


0.1628(49) 


0.3277(71) 


B6 


1.15 


0.1198(52) 


0.6111(91) 


0.627(11) 


0.1313(78) 


0.2066(97) 


B7 


1.175 


0.0660(22) 


0.3593(52) 


0.3659(66) 


0.1083(40) 


0.1055(33) 


B8 


1.18 


0.0565(23) 


0.3085(60) 


0.3199(76) 


0.1108(47) 


0.0927(35) 


B9 


1.185 


0.0430(18) 


0.2292(69) 


0.2277(85) 


0.1090(44) 


0.0664(32) 


BIO 


1.19 


0.0302(16) 


0.1664(81) 


0.165(10) 


0.1083(45) 


0.0506(34) 



TABLE V: PCAC and meson masses from the 24 x 12 3 lattice. 



lattice 


— amo 


am 


aMp S /m 


Mv/Fps 


Mv/Mps 


a 3 {M PS Fps) 2 /m 


BO 


0.95 


0.3931(38) 


5.531(55) 


4.555(81) 


1.03241(85) 


0.618(19) 


Bl 


1 


0.3368(40) 


5.406(67) 


4.637(93) 


1.0376(11) 


0.493(17) 


B2 


1.05 


0.2765(40) 


5.098(77) 


4.75(11) 


1.0428(14) 


0.346(13) 


B3 


1.075 


0.2410(38) 


4.849(78) 


4.85(11) 


1.0421(16) 


0.261(10) 


B4 


1.1 


0.2025(40) 


4.564(93) 


5.03(12) 


1.0418(20) 


0.1810(75) 


B5 


1.125 


0.1604(34) 


4.010(90) 


5.10(14) 


1.0364(30) 


0.1063(51) 


B6 


1.15 


0.1198(52) 


3.12(15) 


4.79(27) 


1.0272(80) 


0.0539(53) 


B7 


1.175 


0.0660(22) 


1.958(74) 


3.38(12) 


1.0183(86) 


0.0229(13) 


B8 


1.18 


0.0565(23) 


1.687(81) 


2.89(13) 


1.036(13) 


0.0207(13) 


B9 


1.185 


0.0430(18) 


1.223(73) 


2.09(11) 


0.993(18) 


0.0145(10) 


BIO 


1.19 


0.0302(16) 


0.918(85) 


1.53(11) 


0.996(39) 


0.01076(98) 



TABLE VI: Mass ratios from the 24 x 12 3 lattice. 
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lattice 


V 


— OWio 




(P) 


T 


A 




CO 


32 x 16 3 


1.15 


5446 


0.665894(44) 


3.32(40) 


0.2227(10) 


3.05(36) 


CI 


32 x 16 3 


1.175 


2192 


0.672235(73) 


2.80(50) 


0.07036(90) 


5.9(1.5) 


C2 


32 x 16 3 


1.18 


4606 


0.673657(49) 


3.46(47) 


0.05167(50) 


6.1(1.1) 


C3 


32 x 16 3 


1.185 


4313 


0.675170(50) 


2.99(39) 


0.03751(38) 


4.66(75) 


C4 


32 x 16 3 


1.19 


5404 


0.676637(44) 


3.29(40) 


0.02474(28) 


7.9(1.5) 



TABLE VII: Bare parameters and average plaquette for the 32 x 16 3 lattice. 



lattice —amo am aMps aMv aFps a 2 Gps 

CO 1.15 0.1175(30) 0.6319(31) 0.6541(43) 0.1196(41) 0.2037(49) 

CI 1.175 0.0678(30) 0.3834(49) 0.4015(61) 0.0919(51) 0.1018(37) 

C2 1.18 0.0549(18) 0.3226(37) 0.3364(46) 0.0860(35) 0.0817(24) 

C3 1.185 0.0420(16) 0.2416(39) 0.2443(50) 0.0784(32) 0.0542(18) 

C4 1.19 0.0308(10) 0.1842(36) 0.1900(43) 0.0806(29) 0.0443(15) 



TABLE VIII: PCAC and meson masses from the 32 x 16 3 lattice. 



lattice 


—amo 


am 


flMp S /m 


Mv/Fps 


My /Mps 


a 3 (Mp S Fps) 2 /m 


CO 


1.15 


0.1175(30) 


3.398(90) 


5.47(18) 


1.0351(35) 


0.0486(26) 


CI 


1.175 


0.0678(30) 


2.17(10) 


4.38(23) 


1.0473(93) 


0.0183(15) 


C2 


1.18 


0.0549(18) 


1.895(70) 


3.91(16) 


1.0426(86) 


0.01404(89) 


C3 


1.185 


0.0420(16) 


1.390(61) 


3.12(13) 


1.011(12) 


0.00854(52) 


C4 


1.19 


0.0308(10) 


1.102(53) 


2.36(10) 


1.031(18) 


0.00715(38) 



TABLE IX: Mass ratios from the 32 x 16 3 lattice. 



lattice 


V 


—amo 


Ntraj 


<^> 


T 


A 




DO 


64 x 24 3 


1.18 


458 


0.673737(46) 


4.0(1.9) 


0.04436(51) 


3.5(1.5) 


Dl 


64 x 24 3 


1.185 


291 


0.675184(59) 


2.3(1.1) 


0.02836(59) 


4.2(2.5) 


D2 


64 x 24 3 


1.19 


349 


0.676649(52) 


1.63(59) 


0.01520(39) 


5.7(3.6) 



TABLE X: Bare parameters and average plaquette for the 64 x 24 3 lattice. 



lattice —amo am aMps aMv aFps a 2 Gps 



DO 


1.18 


0.0562(25) 


0.3433(37) 


0.3597(56) 


0.0620(42) 


0.0661(31) 


Dl 


1.185 


0.0445(36) 


0.2930(73) 


0.325(11) 


0.0629(53) 


0.0580(52) 


D2 


1.19 


0.0330(25) 


0.2184(73) 


0.230(10) 


0.0534(50) 


0.0388(38) 



TABLE XI: PCAC and meson masses from the 64 x 24 3 lattice. 



lattice -amp am aMp S /m My/F PS Mv/M PS a 3 (Mp S F PS ) 2 /m 



DO 


1.18 


0.0562(25) 


2.099(97) 


5.82(37) 


1.0477(98) 


0.00808(93) 


Dl 


1.185 


0.0445(36) 


1.94(19) 


5.20(46) 


1.110(22) 


0.00768(99) 


D2 


1.19 


0.0330(25) 


1.45(13) 


4.34(40) 


1.054(24) 


0.00415(73) 



TABLE XII: Mass ratios from the 64 x 24 3 lattice. 



